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ABSTRACT 


The conversational use of computers makes user -oriented languages desirable. A 
computer program has been developed which enables users to do matrix calculations 
conversationally, without 'knowledge of computers or programming languages. The 
program converses in a brief command language, which provides for input/output, 
matrix arithmetic, forming of the transpose, inverse, determinant, and eigenvalues, 
and exponentiation. Details of the program language and logic are given, with a 
sample user -program conversation. Technical details and FORTRAN program listings 
are also presented. 
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SUMMARY 

Calculations with matrices often require the aid of computers. The ability to do 
such calculations while conversing with the computer frees users from doing any pro- 
gramming to perform everyday calculations, and allows easy checking of results. For 
such an application a user -oriented language, not a computer language, is desirable. 

A computer program has been written which enables users to do matrix computations 
without knowledge of computers or programming languages. Potential users can easily 
learn the command language recognized by the program. The program, called MATAR, 
interacts with users in a conversational manner. It is now running under a time -sharing 
system. 

MATAR provides for simple input/output of matrices, matrix arithmetic, and 
several other operations. These include formation of the transpose, inverse, determi- 
nant, and eigenvalues and exponentiation. 

The limitations of the program are discussed, with some remarks on the logic it 
uses. A complete and brief user’s guide is provided, with a sample user -program 
conversation carried on through a typewriter -like terminal. Technical details and 
FORTRAN program listings are also presented. 


INTRODUCTION 

In the last few years, the increasing capabilities of electronic computers have 
made possible the numerical solutions of complex problems. These problems could 
not practically be solved without computers because of the large amounts of arithmetic 
involved. 

But computers must be told exactly what to do in their own languages. Unfortu- 
nately, this often creates a barrier between the problem solver and his solution. He is 



forced either to master a computer language or to work with a computer programmer who 
translates his problem for the computer. However, there are some problems which 
should not require this; that is, some problems (e.g. , matrix calculations) are just hard 
enough that they require the help of a computer, but not so hard that the problem solver 
should have to translate them into computer language. 

One approach to breaking down this barrier is to develop more capable desk calcu- 
lators which remain easy to use. Desk calculators are not suited to matrix manipula- 
tions, however, because of their storage limitations. 

Another approach is to make computers available for conversational use. By this is 
meant real-time interactions between computers and users such as programmers, busi- 
nessmen, or engineers. The conversational approach is practical especially where time- 
sharing systems are available; under such systems, computers handle more than one 
user essentially as if each user had the computer to himself. 

The conversational approach can be made easy by developing computer programs 
which are strongly user -oriented. For matrix calculations, the user need only know 
how to get matrices in, how to operate on them, and how to get back the answers, all by 
means, for example, of short lines typed at a typewriter -like terminal. This report 
shows a way of doing matrix calculations on a computer by conversing with the computer 
in a brief set of commands. 


MATRIX CALCULATIONS - AN EXAMPLE PROGRAM 
Capabilities of Program 

An unusual feature of MATAR is its power to help users learn the proper usages of 
MATAR. The program can issue specific complaints about input it cannot process. 

For example, if the user tries to use a matrix not previously defined, the program 
tells him that the matrix is unknown and therefore that particular calculation must be 
canceled. The user can promptly see and rectify his errors. 

The program will add, subtract and multiply two matrices, and raise a matrix to a 
positive integer power. In addition, it will form the inverse, determinant, eigenvalues, 
and transpose of user-defined matrices on command. Mixed scalar/matrix arithmetic 
is permitted. The inverse, determinant, and eigenvalues are calculated only for square 
matrices. 

The program stores matrices by names the user chooses and allows him to use 
these names in simple expressions. The user defines a matrix by giving its name and 
a punctuated sequence of numbers for its elements, or by defining it in terms of pre- 
viously defined matrices. The user can quickly display any matrix by asking for it 
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by name. Also, individual elements can be called out or altered. This storage of * 
matrices is limited, but the program allows the user to manage and control use of 
storage by permitting him to delete matrices and to observe what matrices have been 
stored. 

These capabilities make matrix computations quick and simple. Such computations 
have many uses, for example, in the study of linear systems and stability problems. 



(a) First breakdown of user input. 
Figure 1. - General flow of program logic. 
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Limitations of Program 


The program is designed for "medium -size" (10 rows and/or columns) problems. 

It is not intended for solving large systems of equations, for extensive computations, or 
for huge blocks of input or output. For this reason, it is limited in input/ output, 
storage, arithmetic, and reasoning. 

There are two limits on input/ output. First, input/ output is presently by means of 
typewriter -like terminals only, because they are better suited to conversational appli- 
cations than magnetic tape, card readers, or printers. Second, users may type in no 
more than two lines, or about 240 characters. This restriction is easily removed. 

Storage of matrices is limited, both in the size of matrices and in the total number 
of matrix elements. Individual matrices can have as many as 10 rows and/or 10 columns. 
Storage is provided for roughly 1400 elements (e.g. , fourteen 10x10 matrices or fifty-six 
5x5 matrices). 

All arithmetic is limited to six significant figures. Matrix inversion and eigenvalues 
are usually accurate to five figures. The program translates all input matrix elements 
to real, short-precision numbers; only positive integers are allowed in the exponentiation 
operation. 




(b) Utility request processing. 


Figure 1. - Continued. 


(c) Output request processing. 
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Logic limitations allow rapid processing of correct commands, but can occasionally 
cause the user some confusion. For example, the program regards any line with a colon 
on the right side of an equal sign as a definition of a matrix. Figure 1 shows the general 
flow of the logic used to translate commands typed in by the user. This figure is espe- 
cially intended to show why the diagnostic messages are issued. 



(dl Definition request processing. Arithmetic request processing. 

Figure 1. - Continued. 
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(f) Individual element request processing. 
Figure !. ■" Concluded. 


USER'S GUIDE TO MATAR 
Some Details of Usage 

To use the program, users gain access to terminals (or other such devices) and 
begin by executing the call CALL MATAR. Details for doing this will differ (see 
appendix). The following example shows the beginning of a user -program conversation. 
The user commands are (arbitrarily) in small letters; the program responses are in 
capital letters and are indented slightly. Blanks are ignored by the program, except 
that where a number is expected, blanks are interpreted as zero. The user can easily 
include comments because the program ignores lines beginning with an apostrophe. 
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load sample 
£tm sample 

YOU ARE USING A MATRIX ARITHMETIC PROGRAM. 

WANT MORE INFORMATION? 
yes 

I RECOGNIZE THE FOLLOWING COMMAND FORMS: 

(1) NAM El = A11,A12,...,A1M;A21,...,A2M;...;AN1,...,ANM: 

THIS DEFINES MATRIX NAME1, ENTERED ROW BY ROW. 

(2) NAME1 = NAME 2 OP , WHERE OP = \% (INVERSE) D* (DETERMINANT) 

U(TRANSPOSE) PWTO POWER N) OR E%(EI GENVALUES) — 

TH I S LAST DEFINES NAME1 TO BE A MATRIX WHOSE FIRST COLUMN 
IS THE REAL PARTS, AND WHOSE SECOND COLUMN IS THE I MAG 
PARTS, OF THE EIGENVALUES OF MATRIX NAME2. 

(3) NAME1 = NAME 2 OP NAME 3 , WHERE OP = + - * 

(4) NAME1 = ? THIS OUTPUTS MATRIX NAME1 

(5) QUIT, OR STOP THESE CAUSE A RETURN. 

(6) DUMP THIS TELLS YOU ALL MATRIX NAMES AND SIZES. 

(7) ERASE(NAMEl) THIS ERASES MATRIX CALLED NAME1. 

(8) NAM El = NAME 2 OP NAME 3 , WHERE OP = AUDIAGONAL ADD), 

SUDIAGONAL SUBTRACT), OR M|(D I AGONAL MULTIPLY), AND NAME 2 
OR NAME3 MUST BE A SCALAR. 

(9) NAME1(N,M) = XX IF XX IS ? , THIS OUTPUTS THE NAME1 ELEMENT IN ROW N, COLUMN M; 

IF XX IS A NUMBER, IT IS PLACED IN THAT ELEMENT. 

YOU MAY USE TWO LINES FOR INPUT. TO CONTINUE TO A SECOND LINE 
ENTER | SOMEWHERE IN THE FIRST LINE 


DO YOU NEED ADDITIONAL COMMANDS? 
no 

READY 

a=l,0,3.0;0,2,3;0,0,3: 

READY 
ainv=a 12 
READY 
ainv=? 

ROW 1 1.00000 0.0 

ROW 2 0.0 0.500000 

ROW 3 0.0 0.0 

READY 
stop 


- 1.000000 

-0.500000 

0.333333 


CHCIW STOP RETURN 


After MATAR is executing, it asks if information is desired (see preceding example). 
The user replies yes. The program types out information on its usage; when this is 
done, the program types READY. Whenever this READY appears, the program expects 
a new complete command for processing. The example continues by showing the defini- 
tion of a matrix, inversion of it, writing the inverse, and stopping. 

The same problem is done, but with some user’s errors to show typical program 
responses to imperfect commands in the following example: 
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READY 

a=l/ 0, 0, 22,3; 0, 0,3,1*: 

BAD CHARACTER % . ROW 1 , COLUMN 3 

RE-ENTER NUMBER FOR THIS ENTRY. DEFAULT: ZERO. 
3 « 


WARNING: 

READY 

a=? 


ROW 3 IS TOO LONG AND HAS BEEN TRUNCATED. 


ROW 1 1.00000 

0.0 

3.00000 

ROW 2 0.0 

22.0000 

3.00000 

ROW 3 0.0 

READY 
a ( 2 , 2 ) =2 
READY 

aJnv=aa i? 

0.0 

3.00000 

AA NOT FOUND. 

READY 
a i nv=ai % 

READY 
a? n v— i 


DID NOT SEE A ? : 

READY 
a! nv=? 

OR VALID COMMAND OR OPERATOR 

ROW 1 1.00000 

0.0 

-1. 000000 

ROW 2 0.0 

0.500000 

-0.500000 

ROW 3 0.0 

READY 

i dent i ty3«a*al nv 

0.0 

0.333333 

NAME HAS MORE THAN 8 
READY 

i dent 1 ty=a*a i nv 
READY 
ident ? ty=? 

CHARS. BEG 1 

INS WITH IDENTITY 

ROW 1 1.00000 

o 

» 

o 

0.0 

ROW 2 0.0 

1.00000 

0. 596046E-07 

ROW 3 0.0 

READY 
stop 

CHC 1 W STOP RETURN 

0. 0 

1.000000 


LINE REJECTED. 


Command Language 


The program recognizes, and converses in, its own language. This language is 
simple to learn and use. Only this language is recognized, so that knowledge of standard 
computer programming languages is not necessary . 

The user commands are described in the following formats and examples: 

(1) name^ = a^, a 12 , • • *, a i m » a 21» * * *» a 2m’ ■ * •» a nl» * * ’’ a nm : 

This command is used to define an nxm matrix and store it as matrix name j. The 
entries are separated by commas; the rows are separated by semicolons; and the matrix 
must end with a colon. (Note: e-1 is FORTRAN notation for xlO . ) Name^ may con- 
tain one to eight alphanumeric or special characters, with the exception of reserved 
characters, which include + -*%?: = (). 

Example: a - 1, 3.0, 3.2; 1 e-1, 2, 0: 
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This defines matrix a to be 


(2) name^ (i, j) ~ xxx 

This command is used to redefi 
to be the number xxx. 

Example: a(l, 2) = 0 

This redefines matrix a to be 


/l.o 

3.0 

3.2\ 

\0.1 

2.0 

0.0/ 

e the entry in 

the i th row, j th column, 

/l.O 

0.0 

3.2\ 

\o. 1 

2.0 

0.0/ 


(3) name^ = ? 

This command is used to type out the matrix called name^. 

Example: a = ? 

This causes the output 

ROW 1 1.00000 0.0 3.20000 

ROW 2 0.10000 2.00000 0.0 

(4) name.^ (i, j) = ? 

2.1 

This command is used to output the i row, j n column element in matrix name 1 . 
Example: a(2,2)=? 

This command causes the output 2. 00000. 

(5) name^ = nameg op name^ 

where op (short for matrix operator) is + (add), - (subtract), or * (multiply). This 
command is used to perform the arithmetic on the matrices nam^ and name^ and defines 
the result to be matrix name^. K either matrix is lxl, the arithmetic performed is 
scalar arithmetic. 
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Example: b = a + a 

This defines the matrix B to be 


/2.0 

0.0 

6.4 

\0. 2 

4.0 

0.0 


(6) name^ = name 2 op 

where op is i% (inverse), d% (determinant), t% (transpose), p%N (to power N), or e% 
(eigenvalues). The eigenvalue command is used to define matrix name^ to be a two- 
column matrix whose rows contain the real, then the imaginary, parts of the eigenvalues 
of matrix name 2 . 

Example: b = b t% 


This redefines b to be 



As this example shows, previously -defined matrices can be redefined. 

( 7 ) name^ = name 2 op name^ 

where op is a% (add), s% (subtract), or m% (multiply). This command is used to 
perform scalar arithmetic on diagonal entries only. Therefore, either name 2 or name 3 
must be scalar, the other must be square. 

Example: If 



and scalar = 2, c = c m% scalar redefines c to be 
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(8) dump 

This command is used to output all matrix names and dimensions. 

Example (after defining only a and b above): DUMP 
This would cause the output 

YOU HAVE DEFINED 2 MATRICES. 

A HAS 2 ROWS, 3 COLUMNS. 

B HAS 3 ROWS, 2 COLUMNS. 

(9) erase (name^) 

This command is used to erase the matrix name^, releasing the space it occupied for use. 
Example: erase (a) 

This command erases the elements in matrix a. 

(10) stop 

This command stops execution of MATAR. 


PROGRAM MESSAGES 

Whenever communication breaks down and MATAR has trouble with a command, it 
will issue specific complaints. These help the user find errors and learn to use the 
program. There are three basic types of complaints. 

First, no arithmetic is performed (but no previous matrices are lost) if the user’s 
command violates matrix algebra requirements. For example, the program does not 
try to add two matrices with different numbers of rows. Such errors are arithmetic 
errors. 

Second, in trying to fit the input into one of the recognized command formats, the 
program sometimes finds forms which violate built-in assumptions and limitations. 

For example, no name can be longer than eight characters. These violations are called 
user errors. 

These two error types are called fatal, because user commands containing them 
are not carried out. Finally, there are some messages which are for user’s informa- 
tion. They warn him of conditions which are not fatal to the command, but which 
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should be considered. For example, if a user includes nonnumeric characters in a 

number, he must reenter that number, but the rest of the line is not lost. All three 

r 

types of program messages are listed in table I. 


TABLE I. - MESSAGES ISSUED BY PROGRAM 


User errors 


1 NAME HAS MORE THAN 8 CHARACTERS. BEGINS WITH xxxxxxxx. 

2 RESERVED CHARACTER USED IN NAME, LINE REJECTED. 

3 xxxxxxxx NOT FOUND. 

4 BLANK(S) FOUND WHERE NAME, COMMAND, OR REPLY REQUIRED, RE-ENTER. 

5 DID NOT SEE A ? : OR VALID COMMAND OR OPERATOR, LINE REJECTED. 

6 MORE THAN 10 COLUMNS, LINE REJECTED. 

7 MORE THAN 10 ROWS, LINE REJECTED. 

8 SUBSCRIPT(S) OUT OF RANGE, xxxxxxxx HAS m ROWS, n COLUMNS. 

9 BAD CHARACTER x IN EXPONENT, LINE REJECTED. 

10 BAD CHARACTER x IN NUMBER, LINE REJECTED. 

11 BAD CHARACTER x IN SUBSCRIPT, LINE REJECTED. 

12 NO COMMA SEPARATING SUBSCRIPTS, LINE REJECTED. 

Arithmetic errors 


13 MATRICES HAVE DIFFERENT NO. OF ROWS. OPERATION IGNORED. 

14 MATRICES HAVE DIFFERENT NO. OF COLUMNS. OPERATION IGNORED. 

15 OPERATION NOT DEFINED FOR NONSQUARE MATRIX. 

16 FIRST NO. COLUMNS NOT EQUAL SECOND NO. ROWS. 

MATRIX MULTIPLY REQUEST IGNORED. 

17 SINGULAR MATRIX. INVERSE NOT PERFORMED. 

18 OPERATOR NOT APPLICABLE FOR MATRICES. 


Warnings 

19 MATRIX ILL-CONDITIONED. INVERSE MAY BE INACCURATE. 

20 * * * WARNING * * * ROOM FOR ROUGHLY n MORE ENTRIES. 

21 *****WARNING***** STORAGE SPACE FOR MATRICES EXCEEDED BY n, 
RESULTS UNPREDICTABLE. 

22 BAD CHARACTER x. ROW n, COLUMN m. 

RE-ENTER NUMBER FOR THIS ENTRY. DEFAULT: ZERO. 

23 WARNING: ROW n IS TOO LONG AND HAS BEEN TRUNCATED. 

24 WARNING: ROW n IS TOO SHORT AND HAS BEEN PADDED WITH ZEROS. 

25 **WARNING** EIGEN FAILED TO CONVERGE ON AN EIGENVALUE. 
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CONCLUSIONS 


Using computers conversationally, that is, in a real-time interactive way, can be 
easy for nonprofessionals, as well as fruitful for computer professionals. Matrix cal- 
culations have been presented as an example of this conversational use. A brief but 
powerful set of commands to do matrix calculations has been developed. A FORTRAN 
program has been written which enables the computer to converse with users (for 
example, through typewriter -like terminals) and interpret these commands. This 
program is running successfully under a time-sharing computer system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 27, 1968, 

125-23-02-15-22. 
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APPENDIX - NOTES FOR PROGRAMMERS 


Program Details 

The program was developed and run under the Time-Sharing System (TSS) of an IBM 
360/ model 67 (ref. 1). The entire program is roughly 9200 four -byte words (1 byte = 8 
bits) in length. Under TSS, each user’s virtual storage is divided into pages, at 4096 
bytes per page. The number of these pages that the program requires is a strong 
function of how much the various subprograms have been combined at compile time and 
by the link-editor. 

Source decks for the program are in TSS FORTRAN (ref. 2). The program calls two 
shift functions standard to IBM 360 - ISLL (NN, IWORD), which shifts the integer short- 
precision word IWORD NN bits to the left- and ISRL(NN, IWORD), which shifts IWORD NN 
bits to the right. (See note below. ) The program is executed under TSS by compiling the 
source decks, and then LOADing and RUNning the object decks into the user’s virtual 
memory. 

Note: H shift functions above are not readily available, the following FORTRAN coding 
will serve: 


C FORTRAN SHIFTS FOR MATAR 

FUNCTION ISLLCNN, IWRD) 

C MASKR GOOD ONLY FOR 24 BITS TO RIGHT 

INTEGER MASKR/Z00Q00080/, MASK/Z80000000/ 
IWORD = IWRD 
DO 1 I =1, NN 
IWORD = I W0RD*2 

1 CONTINUE 
GO TO 3 

ENTRY ISRLCNN, IWRD) 

IWORD = IWRD 

C SET SIGN BIT TO ZERO 

I F ( I WORD. LT. 0) IWORD = IWORD + MASK 
DO 2 1=1, NN 
IWORD = I WORD/2 

2 CONTINUE 

C RESTORE BIT IF NEEDED 

IF< IWRD.LT.0) IWORD = IWORD + MASKR 

3 ISLL = IWORD 
RETURN 

END 
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Numerical Methods 


Matrix arithmetic is done by real short-precision arithmetic with operations t>n 
32 -bit words yielding six significant figures. 1%© exponentiation command simply causes 
repeated multiplication. 

The inverse, determinant, and eigenvalues are calculated only for square matrices. 
The inverse and determinant are formed by using Gauss -Jordan elimination (ref. 3). The 
eigenvalues are calculated by transforming the matrix using Danilevski's method (ref. 4). 
These calculations are performed by means of standard subprograms developed at this 
Center. 


Listing of FORTRAN Program 

The program is entered by running any program which makes the call CALL MATAR. 
All matrix information and elements are in COMMON/ MATE / . A complete listing 
follows. 


INITIALIZES* TALUS A LITTLE 
SUBROUTINE MATAR 
DIMENSION INPI1G) ,NM(2) 

COMMON/ MATE/' N.UM , 1 L AST , MATR I X ( 1500) 

DATA lY,jN,INEW/» YES',* NO‘,1/ 

IFUNEW.N6.1 > GO TO 3 

I NEW ® 2 

NUM = 0 

I LAST = 0 

DO 5 1*1,1500 

5 MATRIX ID = 0 
PRINT 101 

GO TO 2 

1 PRINT 102 

2 READ 201 » IMP 

CALL MM PACK ( INP, 1 , 10,NM, ,TRUE . ,E1 ,61,62) 

IF<MM(2').E0,IHt GO TO 3 

If (NM<2>.NE.IY) GO TO 1 

PRINT 103 

PRINT 106 

PRINT 104 

GO TO 7 

6 PRINT 102 

7 READ 201, IMP 

CALL UNPACK U HP , 1 , 10 ,W , , TRUE » , 6.6 , 6.6 » 67 ) 
IFINM2) .EQ.1N) GO TO 3 
I F I NM(2 ) »NE. I Y > GO TO 6 
PRINT 105 

3 CALL READER 
RETURN 
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101 FORMAT { * YOU ARE USING A MATRIX ARITHMETIC PROGRAM.*,/, 

4 I WANT MORE INFORMATION* 1 * ) 

102 FORMAT < * TYPE YES OR NO,*) 

103 FORMAT { * J RECOGNIZE THE FOLLOWING COMMAND FORMS,*,/* 

1 ‘ (1) NAME 1 « AH*A12*,..*AlMjAPlt...*A2M j...j ANl*,,,*ANM r *,10X, 

2 / *5X, ' THIS DEFINES MATRIX NAME! * ENTERED ROW BY ROW.',/* 

3 * (2) NAME! = NAME2 OP * WHERE OP * 1% (INVERSE) 0% (DETERMINANT 

4)>,/,5X t * W TRANSPOSE) PINT TQ POWER N) OR E&I 6 1 GENVALUES ) — * * / » 
CSX,* THIS LAST DEFINES N&Mfil TO BE A MATRIX WHOSE FIRST COLUMN*,/ 
D, 5X , ’ IS THE REAL PARTS* AND WHOSE SECOND COLUMN IS THE JMAG * , / , 
E5X,» PARTS, OF THE EIGENVALUES OF MATRIX NAME 2* * , / * 

5 * (3) NAME! = NAME 2 OP NAME 3 * WHERE OP ~ + ~ * **/* 

ft * (4) NAME1 * 7 THIS OUTPUTS MATRIX NAME!',/, 

T* (5) QUIT, OR STOP THESE CAUSE A RETURN,*,/, 

ft » (6) DUMP THIS TELLS YOU All MATRIX NAMES AND SIZES,*,/, 

9 * (?) ERASE (NAMED THIS ERASES MATRIX CALLED NAME 1 . * ) 

104 FORMAT (» DO YOU NEED ADDITIONAL COMMANDS?*) 

105 FORMAT (' CONTACT BRUCE C.ANRIGHT, PAX 6263 OR 8224,*) 

106 FORMAT (* (ft) NAME! * NAME2 OR NAME3 , WHERE OP ~ A$>( DIAGONAL ADD) 
1,* ,/,5X,‘S^iDIA60NAL SUBTRACT), OR MSS (DIAGONAL MULTIPLY), AND NAME 
22 * */ ,5X , * OR NAME 3 MUST BE A SCALAR. »,/,» (9) NAME 1 ( N*M) » XX 

3 IF XX IS 7 * THIS OUTPUTS THE NAME 1 ELEMENT JN ROW N, COLUMN M ? »,/ 
3, 5X, * IF XX IS A NUMBER, IT IS PLACED IN THAT ELEMENT.*,/,/, 

A » YOU MAY USE TWO LINES FOR INPUT, TO CONTINUE TO A SECOND LINE* 
B */,* ENTER { SOMEWHERE IN THE FIRST LINE*,/,/ ) 

201 FORMAT ( 10. A 1 ) 

END 


C 


SUBROUTINE READER 

: # ft#.#.## # ajtsje # & jjc # ^ if* £#### # # ## # Sjc# % & 3p£ # # % %% & ^ & & # # & #s 


READER READS THE INPUT AND DECIDES WHAT TO CALL 
THE 1,2 ENTRY IS IN THE FIRST ROW, SECOND COLUMN 
OUTLINE OF MAJOR FLOW FOLLOWS 

FOR COMMAND TYPE (li RCADER— CHANGE — CNVTR 
FOR COMMAND TYPE (2), (4) READER — REPAIR 
FOR COMMAND TYPE READER— -DOT PUT 

FOR COMMAND TYPE READER — ARITMM — AR I 7HMETICSTUFF 

FOR COMMAND TYPE (8), (9), (ID) READER — COMANQ — DUMPY 

NMPACK PACKS A*1 LEFT-ADJUSTED INTO 2A4 RIGHT-ADJUSTED 
NMDEF DEFINES NEW NAME 
FINDNM LOOKS UP OLD MATRIX 
WIPE ERASES OLD MATRIX 




IS 
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LOGICAL OLD 

DIME N SION N.M < 2 ) , N M 1 ( 2 ) , I NPUT (220) 

COMMOM/MATE/MOM, I LAST , HATH I X { 1500) 

C 0 M M ON / PASS./ K L P , K R P , K , K C 0 L 

DATA I I , ID, IT, IP, IE / * I 1 , 1 D 1 , « T * t • % «, 1 P 1 / 

DATA I.M, IA, IS, ICON, IEO, IQM, IBLK/ 1 * 1 , • + • , '- 8 , 1 f 1 , 8 = 1 * 8 V » 8 •/ 

DATA I C OM , I S EM , I COL , 1 L P , I R P , I V / 1 , 8 , 8 •* 9 1 : 8 , 1 < 8 , 1 ) 8 1 8 £ 8 / 

DATA I A 2 , I S 2 , I M 2 , I A P / ’AS^SS 8 ^, 1H* / 

50 IFULAST.LT.1400) GO TO 47 
NMA.X = I LAST - 1399 
PRINT 98, NMAX 

98 FORMAT! 1 ❖ ** * * WARN IMG ** * * * STORAGE SPACE FOR MATRICES EXCEEDED*,/ 

*, 8 BY * , I 3 , 1 , RESULTS UNPREDICTABLE . * ) 

GO TO 46 

47 IF( ILAST.LT. 1300) GO TO 46 
MM AX - 1399 - I LAST 
PRINT 99, NM AX 

99 FORMAT! » ***WAR-NING*** ROOM FOR ROUGHLY *,I4,* MORE ELEMENTS. * I 
46 PRINT 100 

100 FORMAT! « READY*) 

DO 45 1=1,220 
45 INPUT! I) = IBLK 

READ 201, ( INPUT! I) ,1=1,110) 

DO 51 1=1,110 
K = I 

IF! INPUT! K) .EQ. ICON) GO TO 52 

51 CONTINUE 

201 FORMAT ( 120A1 ) 

GO TO 53 

C LINE CONTINUED 

52 READ 201, (Z NPUT ( I ) , I =1 1 1 , 220 ) 

INPUT (K) = IBLK 

LOOK FOR EQUALS SIGN 

IF APOSTROPHE, IGNORE WHOLE LINE, ASSUME COMMENT 

53 DO 54 1=1,220 
KEQ = I 

IF ( INPUT ( KEQ ) . EQ • I AP ) GO TO 46 
IF! INPUT { KEQ ) . EO. I EQ ) GO TO 55 

54 CONTINUE 

NO EQUALS SIGN IS IT COMMAND 
CALL C OMAND ( I N P U T , NM , £60 ,£5 , £44 , £46 ) 

GO TO 46 
44 RETURN 

LOOK FOR LHS PARENS 

55 DO 551 1=1, KEQ 
KLP = I 

IF! INPUT(KLP) .EQ. LLP) GO TO 552 
351 CONTINUE 
GO TO 560 

552 00 553 I =KLP ,KEQ 
KRP = I 

IF! INPUT(KRP) .EQ.IRP) GO TO 554 

553 CONTINUE 
GO TO 49 

554 DO 555 I=KLP,KRP 
KCOM = I 

IF! INPUTIKCOM) .EQ. IC.OM) GO TO 556 
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VBbVjTNTINUE 
PR l NT 101 

101 ? FORMAT ( » NO COMMA SEPARATING SUBSCRIPTS* LINE REJECTED.') 

* GO TO 46 

556 CALL NMPACK ( I NPUT , l ,KLP-l ,NM t .F ALSE. , £ 56 , £49, £46 ) 

CALL FINDNM(NM, INM,£60) 

CALL REPAIR (INPUT, I MM ,KLP*KEQ*KRP* NM *646 ) 

GO TO 46 

C ASSUME REGULAR MATRIX NAME LHS 

C PACK THE LEFT HAND NAME 

560 INPUT ( KEQ ) = I BLK 

CALL NMPACK ( I NPUT , 1 , KEO , NM , .FALSE. ,556* 649*546) 

GO TO 57 

56 PRINT 106, NM 

106 FORMAT ( 1 NAME HAS MORE THAN 8 CHARS. BEGINS WITH *,2A4) 

GO TO 46 

49 PRINT 105 

105 FORMAT ( ' RESERVED CHARACTER USED IN NAME, LINE REJECTED.') 
GO TO 46 

C LOOK FOR QUESTION MARK 

57 DO 58 I =K EQ , 220 

IF( INPUT ( I ) .EQ. IOM) GO TO 59 

58 CONTINUE 

GO TO 90 

C ASSUME NAME = ? 

59 CALL F I NDNM ( NM * I NM * £60 ) 

C NAME FOUND 

NROWS = MATRIX ( INM+3) 

NCOLUM = MATRIX! INM+4) 

CALL OUTPUT (MATR IX( I NM+5 ), NROWS , NCOLUM ) 

GO TO 46 

C NAME NOT FOUND 

60 PRINT 107, NM 

107 FORMAT ( 5X,2A4,' NOT FOUND.') 

GO TO 46 

61 PRINT 107, NM1 
GO TO 46 

62 PRINT 106, NM 1 
GO TO 46 

65 DO 66 I=KEQ*220 

C LOOK FOR VALID OPERATORS 

K = I 

IF( INPUT (K) .EO.IM) GO TO 71 
IF( INPUT (K ) .EQ. IA) GO TO 72 
I F ( I NPUT (K).EO.IS) GO TO 73 

66 CONTINUE 

C NO * + OR - ENCOUNTERED. LOOK FOR INV OR DET. 

C LOOK FOR I % OSg T% 

DO 67 I =KEO * 2 1 8 
K = I 

I F ( INPUT ( K ) .EQ. IP) GO TO 68 

67 CONTINUE 
GO TO 5 
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6R K = K - 1 

HEMP = INPUT (K) 

IFMTEMP *B0. II) GO TO 74 
I F ( I TEMP .EQ. ID) GO TO 75 
I F ( I TEMP * EO * IT) GO TO 69 
IF (I TEMP * EO * IV) GO TO 70 
IF l I TEMP. EO. IE) GO TO 81 
I F ( I TEMP • EQ * I A 2 ) GO TO 10 
I F ( I TEMP. EO * I S2 ) GO TO 11 
IF ( ITEMP.EQ.IM2) GO TO 12 
C NO KEY THINGS FOUND 

5 PRINT 109 

109 FORMAT ( 1 DID NOT SEE A ? * OR VALID COMMAND OR OPERATOR. LINE 

REJECTED. 1 ) 

GO TO 46 

C OPERATOR CODES * 1 + 2 - 3 INV 4 DET 5 

69 KODE = 6 
GO TO 76 

70 KODE a 8 
GO TO 76 

71 KODE = 1 
GO TO 76 

72 KODE = 2 
GO TO 76 

73 KODE =‘ 3 
GO TO 76 

74 KODE = 4 
GO TO 76 

75 KODE = 5 
GO TO 76 

TO KODE = 9 
GO TO 76 

11 KODE = 10 
GO TO 76 

12 KODE =11 
GO TO 76 

81 KODE = 7 
KJ = K + 1 
INPUT ( K J ) = ILP 
INPUT (220) = IRP 

CALL CNVTR ( I NPUT , K J ,220 , 2 , ARC , I NM2 1 1 CHAR , £.76 , S83 * £84 ) 

GO TO 76 

83 I CHAR = ICON 

84 PRINT 121 t I CHAR 

121 FORMAT ( 1 BAD CHARACTER ',A4,' IN EXPONENT, LINE REJECTED.') 

GO TO 46 

C CHECK TO SEE IF LEFT-HAND NAME IS NEW 

76 OLD = .FALSE. 

CALL F I'NONM { NM , INM ,879 ) 

OLD- = .TRUE. 

79 J J = KEO + 1 

C RESOLVE R.IGHT.HA.MD NAME (S'). 

I NPUT ( K ) = IBLK 

CALL NMPACK ( INPUT , J J , K , NM 1 , .TRUE . ,862 , 849,846) 

CALL FINQNMi NM1 , INM1 ,861 ) 

NA ?s. MATRIX ( INM1+.3) 

MA = MA TR I X ( INM'1+4) 

I F ( KODE# GT .3 .AND. KODE .LT .9 ) GO TO 77 
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JJ = K+l 

I F { KODE . GT • 3 ) JJ = JJ + 1 

CALL NMPACKUNP'UT, JJt220 f NMl f /TRUE, * £62 T E49, S46 ) 
CALL FINDNM(NM1,INM2,£61> 

MB = MATRIX ( I NM2+3 ) 

MB = MATRIX ( I MM2 +4 } 

C PUT AWAY NEW NAME 

77 CALL NMDEF(NM) 

C NAMES RESOLVED, READY TO DO ARITHMETIC 

IN - I LAST - 1 

C ALL AR I TH (KODE * I N , I MM 1 , I MM2 , NA , MA , MB , MB , NC , MC , £80 ) 
MATRIX! ILAST+2) = NC 
MATRIX { I L AST + 3 ) = MC 
MATRIX* ILAST+1 ) = NC*MC + 5 
I LAST = ILAST + NC*MC + 3 
I F ( .NOT ® OLD ) GO TO 50 
C ERASE OLD NAME IN MIDDLE 

CALL WIPE( INM) 

GO TO 50 

C BAD MATCHUP ROWS AND COLUMNS GO BACK 

80 MATRIX* ILAST) = IBLK 
MATR I X ( I LA'S T-l } = IBLK 
ILAST = ILAST - 2 
NUM = NUM - I 
GO TO 46 

C LOOK FOR COLON IF SO DEFINE MATRIX 

90 DO 91 I =KEQ ,190 
KCOL = I 

I F ( INPUT ( I ) • E Q * ICOL) GO TO 92 

91 CONTINUE 

C NO COLON FOUND ASSUME A-RITH 

GO TO 65 

92 INPUT(KEQ) = ILP 
KLP ~ K EO 

KRP = KCOL + 1 
INPUT(KRP) = IRP 

95 NC DM = 0 

CALL F I NDNM { NM , I NM , £96 ) 

C ALREADY DEFINED 

CALL W I PE ( INM) 

96 CALL NMDEF(NM) 

C COUNT THE FIRST COMMAS , FIND A SEMICOLON AT K * 

DO 1 I =KL P , KCOL 
K = I 

IF * INPUT (K) * EQ • I COM) NCOM = NCOM + 1 
I F ( I N P LJT { K ) # EQ • I S EM *0R. INPUT ( K ) . EQ • I COL ) GO TO 2 

1 CONTINUE 

2 NCOLUM = NCOM + 1 
NSEMI = 0 

I F ( INPUT ( K ) • EO • I SEM ) NSEMI = 1 
I F ( NCOLUM « L E • 10) GO TO 21 
PRINT 102 

102 FORMAT ( 1 MORE THAN 10 COLUMNS • LINE REJECTED*') 

GO TO 80 

21 INPUT(K) = IRP 
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C COUNT THE SEMICOLONS T COLON AT J . 

1)1) 3 I =K t KCUL 
J = I 

I F ( I N P U T ( J ) .EO. I S EM ) NS EM I = NS EM I + l 
3 CONTINUE 
A NROWS = NSEM1 + 1 

I F ( NROWS .LE . 10) GO TO 41 
PRINT 103 

103 FORMAT ( * MORE THAN 10 ROWS. LINE REJECTED.*) 

GO TO 80 

41 MATRIX! ILAST+1 ) = NROWS*NCOLUM + 5 
MATRIX! I LAST +2 ) - NROWS 
MATRIX! ILAST+3) = NCOLUM 
I LAST = I LAST + 3 

CALL CHANGE ( INPUT .MATRIX! ILAST+1) , NROWS .NCOLUM, &97 ) 
I LAST = I LAST + NROWS*NCOLUM 
GO TO 50 

C GO BACK , CONVERSION FAILED . 

97 I LAST = I LAST - 3 
GO TO 80 
END 



SUBROUTINE NMP ACK ( I NPUT , K1 , J 1 » NM .OLD , * 


: # 3' t c sjc s;? ^ 3fc sjs # ❖ ❖ ^ 


ISLL A NO ISRL ARE ASSEMBLY LANGUAGE ROUTINES IN OUR SYSTEM 
ISLLINN, I WORD) SHIFTS 1*4 IWORD NN BITS TO LEFT 
ISRL ( NNf I WORD ) SHIFTS 1*4 IWORD MM BITS TO RIGHT 


< : 


c # >\z # jJc sjc sjc-sjc sjc # s£ # sjc >;? y/: sjc ;; 


TAKES THE JUNK IN IMPUT(K) TO INPUT! J ) AND 

PACKS IT INTO A NAME IF IT LESS THAN 9 NON-BLANK CHARACTERS 
RETURN 1 MORE THAN 8 CHARS 
RETURN 2 OPERATOR IN NAME 
RETURN 3 BLANKS FOUND 
LOGICAL OLD 

DIMENSION NM ( 2 ) .INPUT ( 1 ) 

DATA I BL , I PL , I MI , I AS , I PER , I OM , ICOL, I EG, ILP, IRP 

/ • , • % • t » , ' = • ' / 

j = ji 
K = K 1 
NM ( 1 ) = I BL 
NM (2 ) = I BL 
KTR = 0 

IF (OLD) GO TO 3 
DO 1 I— K , J 
IN = INPUT! I ) 

IF! IN. EO . I PL .OR. IN.EQ.IMI .OR. IN.EO.IAS .OR. I N. EO . I PER ) GO TO 6 
IF! IN. EO. ILP .OR. IN. EO. IRP) GO TO 6 

IF! IN.F0.IOM .UR. IN. EO. ICOL .OR, IN.EO.IEO) GO TO 6 
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1 continue 

3 no 4 I=K,J 

I TEMP = INPUT (I) 

I F { I TEMP. £0. IBL) GO TO 4 
KTR = KTR + 1 
IF(KTR.GT.B) GO TO 5 

NM ( 1 ) = ISLL! 8,NM( 1 ) ) + I SRL ( 24, NM ( 2 ) ) 

NM ( 2 ) = I SLL ( 8 , NM ( 2 ) ) + I SRL ( 24 , 1 TEMP ) 

4 CONTINUE 

IF(KTR.EQ.O) GO TO 7 
RETURN 

5 RETURN 1 

6 RETURN 2 

7 PRINT 11 

11 FORMAT ( 1 BLANK (5) FOUND WHERE NAME, COMMAND, OR REPLY REQUIRED, R 
’•-E-ENTER . * ) 

RETURN 3 
END 


DEFINE NAME AND PUT IT IN BIG ARRAY 
SUBROUTINE NMDEF(NM) 

DIMENSION NM ( 2 ) , INPUT ( 1 ) 

COMMON/MATE/ NUM, I LAST , MATR I X (1500 ) 

NUM = NUM + 1 
I LAST = ILAST + 2 
MATRIX! ILAST) = NM ( 2 ) 

MATRIX! I LAST-1 ) = NM ( 1 ) 

RETURN 

END 


SEARCH FOR NAME 

RETURN 1 NM NOT FOUND 

RETURN NM ( 1 ) AT MATRIX! I NM ) 

SUBROUTINE F I NDNM ( NM » I NM , # ) 

DIMENSION NM ( 2 ) 

COMMON/MATE/ NUM, I LAST ,MATR I X ( 1500 ) 

J = 1 

DO 10 1=1, NUM 

IF ( NM ( 1 ) .EQ.MATR I X ( J ) .AND. NM ( 2 ) . EG . MATR I X ( J+ 1 ) ) GO TO 15 
J = J + MATRI X ( J + 2 ) 

10 CONTINUE 
RETURN 1 
15 INM = J 
RETURN 
END 
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ERASES MATRIX NM, MOVES EVERYTHING BELOW UP 
SUBROUTINE WIPE! I I) 

COMMON/MATE/ NUM , I L AST , MA TR I X < 1500 ) 

INM = II 

I DEL = MATRIX! I NM+2 ) 

1ST = INM + I DEL 
I END = 1ST - 1 
DO 2 I =INM» I END 
2 MATRIX! I) = 0 
DO 1 I = I ST , 1 LAST 
J = I - I DEL 
MATRIX! J) = MATRIX! I ) 

1 CONTINUE 

I LAST = I LAST - I DEL 
NUM = NUM - 1 
RETURN 
END 


THIS IS THE ONLY WAY TO GET STUFF OUT 
SUBROUT I N E OUTPUT ( ARRAY » NR , NO 
DIMENSION ARRAY (NR, NO 
MROWS = NR 
NCOLUM = NC 

I F ( NROWS .EO. I) GO TO 15 
DO 10 1=1, NROWS 

10 PRINT 101, I , (ARRAY! I, J) ,J=1, NCOLUM) 
RETURN 

15 PRINT 102, ( ARRAY ( 1 ,J ) ,J=1 , NCOLUM) 

RETURN 

101 FORMAT ( 1 ROW' , 1 3 , 5G 14 . 6 , / , 5X , 5G 14. 6 ) 

102 FORMAT ( 5X » 5G14- .6 , / » 5X » 5G14.6 ) 

END 


C 

C 


SUBROUTINE REPAIR ( INPUT ,INM ,ILP , I EG ,KRP ,NM ,* ) 
PRINT OUT OR REPLACE INDIVIDUAL ENTRIES 
RETURN 1 BAD SUBSCRIPT OR NUMBER 
DIMENSION I NPUT ( 1 ) , NM ( 2 ) 

COMMON/MATE/ NUM, I L AS T , MATR I X ( 1500) 

DATA IRP, ILP, I OM , I COM / ' / 

IPT = INM 

NROW = MATRIX! IPT+3) 

NCOL - MATRIX ( IPT+4) 
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C CONVERT THE SUBSCRIPTS 

KKK = KLP 
DO 2 1=1,2 

CALL CNVTR ( I NPUT , KKK ,KRP , 2 , ARG , I ARG , I CHAR , S 1 5 , S 1 , £ 15) 

ICOL = I ARG 
GO TO 3 

1 IROW = I ARG 

2 CONTINUE 

3 IFUROW.LE.O •OR. ICOL.LE.O •OR. IROW.GT.NRQW .OR. ICQL.GT. 

* NCOL) GO TO 14 

C SUBSCRIPTS OK 

C GET POINTER TO PROPER ENTRY 

I PT2 = NROW*{ ICOL- 1 ) + IROW + IPT+U 
C IS THERE A QUESTION MARK 

DO 5 I = I EQ, 190 

I F ( INPUT ( I ) .EQ. IOM) GO TO 9 

5 CONTINUE 

...NO, ASSUME NUMBER 
INPUT (190) = IRP 
I NPUT ( I EQ ) = ILP 

CALL CNVTR ( INPUT , I EQ , 190 , 1 , I ANS , I ARG , I CHAR , L6 , £7 , SB ) 

6 MATRIX! I PT2 ) = IANS 
GO TO 10 

7 I CHAR = ICON 

8 PRINT 103, ICHAR 

103 FORMAT ( * BAD CHARACTER ' , A4 , ' IN A NUMBER, LINE REJECTED.*) 
GO TO 20 

...YES, WRITE IT OUT 

9 CALL OUTPUT (MATRIX! I PT2 ) , 1 , 1 ) 

10 RETURN 

SUBSCRIPTS OUT OF BOUNDS 

14 PR I NT 101, N M , N R 0 W , N'C 0 L 

101 FORMAT ( 1 SUBSCRIPTS) OUT OF RANGE, *,2A4,> HAS ',12, 

* « ROWS, ',12,* COLUMNS.') 

GO TO 20 

BAD SUBSCRIPT 

15 PRINT 102, ICHAR 

102 FORMAT ( ' BAD CHARACTER *,A4,' IN SUBSCRIPT, LINE REJECTED.*) 
20 RETURN 1 

END 


SUBROUTINE CHANGE { INPUT, ARRAY , NROWS , NCOL AM , * ) 

CONVERTS THE INPUT STRING TO R*4 NUMBERS, PUTS INTO /MATE/ 

DIMENSION INPUT (1) , INTEMP ( 90 ) , ARRAY ( NROWS , NCOL AM ) 

COMMON/PASS/ KLP,KRP,K,J 

DATA ILP, IRP, ISEM, ICOM, ICOL / '( T 1 , ' , 1 t ' / 

DATA I NTEMP ( 1 ) , INTEMP ( 90 ) / * <»,')*/ 

NCOLIJM = NCOL AM 
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C CONVERT THE FIRST ROW 

KKK = KIP 
DO 5 1=1, NCOLUM 

CALL CNVTR ! I NPUT , KKK , KRP , 1 , ANS , GARB , ICHAR , £20, £ 10, £4 ) 

C ROW DONE 

ARRAY! 1,1) = A NS 
GO TO 6 

C BAD CHARACTER , RE-ENTER 

1 ICHAR = ICON 

4 PRINT 103, ICHAR, I 
PRINT 105 

READ 201, UNTEMP! I I I ), I I 1=2 ,89) 

L = 1 

CALL- CNVTR! I NTEMP , L , 90 , 1 , ANS ,GARB , ICHAR , £20 , £ 1 , £4) 

C SKIP ON TO NEXT COLUMN (THANKS DUE TO EVA) 

DO 2 I 2 = KKK , KR P 
ITT = INPUT! 12) 

IF! I TT • EO . ICOM) GO TO 3 

I F ( ITT • EO . I SEM .OR. ITT.EQ.ICOL) GO TO 6 

2 CONTINUE 

3 KKK = 12 

C NUMBER DONE, ROW NOT YET DONE 

10 ARRAY (1,1) = ANS 

IF! I • EO . NCOLUM ) GO TO 6 

5 CONTINUE 

20 PRINT 102 

102 FORMAT! 1 NO RIGHT PAREN. SOMETHING REALLY WRONG.’) 

RETURN 1 

103 FORMAT (« BAD CHARACTER ' ,A4,«. ROW . 1 , COLUMN * , 13 ) 

105 FORMAT! 1 RE-ENTER NUMBER FOR THIS ENTRY. DEFAULTS ZERO. 1 ) 

201 FORMAT { 8 8A1) 

C DO ROW SUB KTR 

6 KTR = 2 

IF(NROWS.EO, 1 ) RETURN 

7 INPUT ! K ) = ILP 
NCOM = 0 

DO 8 I =K , J 
K2 = I 

IF! INPUT (K2) * EQ # I SEM.OR • INPUT (K2 ) .EO.ICOL) GO TO 9 
IF! INPUT (K2) .60. I COM) NCOM = NCOM + 1 

8 CONTINUE 

C FOUND NEXT SEMICOLON , MOVE RIGHT TO K2, SAVE LEFT K 

9 NMAX = NCOM + I 

IF (NCOLUM-NMAX ) 23,11,21 

21 PRINT 107, KTR 

107 FORMAT! 1 WARNINGS ROW 1 , 13,* IS TOO SHORT AND HAS BEEN PADDED WITH 
* ZEROS.*) 

C PAD IT 

NMAX = NMAX + 1 
DO 22 I = NM A X , NC OL UM 

22 ARRAY ! K TR , I ) = 0.0 
GO TO 11 

23 PRINT 104, KTR 

104 FORMAT!' W ARN I NG * ROW 1 , I 3 , 1 IS TOO LONG AND HAS BEEN TRUNCATED.') 

C CONVERT THE ROW 

11 INPUT ! K2 ) “ IRP 
KKK = K 

DO 15 1=1 , NCOLUM 

CALL ' CNVTR! I NPUT , KKK , KRP , 1 , ANS , GARB , ICHAR , £20 , £12 , £ 14 ) 
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C ROW DONE 

ARRAY ( KTR * I ) = ANS 
GO TO 16 
24 ICHAR = ICOM 

14 PRINT 106 , ICHAR, KTR, I 

106 FORMAT ( * BAD CHARACTER • ,A4,' . ROW*, 13,' COLUMN', 13 ) 

PRINT 105 

READ 201, < INTEMP! 1 1 1 ), I I 1=2,89) 

L = 1 

CALL CNVTR ( IMTEMP,L, 90,1, ANS, GARB, ICHAR, £20, £24,614) 

C SKIP ON TO NEXT COLUMN 

DO 17 I 2=KKK , KRP 
ITT .= INPUT! 12 ) 

I F ( ITT. EQ. ICOM) GO TO 18 

I F ( ITT. EO. I S EM .OR. ITT.EO.ICOL) GO TO 16 

17 CONTINUE 

18 KKK = 12 

C NUMBER DONE, ROW NOT YET DONE 

12 ARRAY ( KTR , I ) = ANS 

15 CONTINUE 

16 KTR = KTR + 1 
IF(KTR.GT . NROWS ) RETURN 
K = K2 

GO TO 7 
END 


C 


SUBROUTINE CNVTR ( INPUT , KTR, MAX, TYPE, ARGl ,ARG2, CHAR, *,*,*) 


Q >Jc if * ** $ ncsjc 


* if if if if i^ if if if if if if if if if if * * * >5= if if * * 'if. if if if if if if if # * if if if 5jc if * if if if ^ if if if if * if ** if ; 


SUBROUTINE TO CONVERT A TYPE INPUT TO FLOATING POINT NUMBERS 

INPUT ARRAY OF CHARACTERS • ONE PER WORD LEFT AOJ. 

KTR INDEX OF FIRST WORD OF ARRAY TO BE CONSIDERED 

MAX DIMENSION OF INPUT 

TYPE =1 REAL*4, =2 INTEGER*^ 

ARGl REAL-4 ANS 

ARG2 INTEGER-4 ANS 


CHAR 
RETURN 
RETURN 1 
RETURN 2 
RETURN 3 


BAD CHARACTER IF FOUND 
NEXT CHARACTER ) RESULT GOOD 
NO , OR ) FOUND RESULT ZERO 
NEXT CHARACTER * RESULT GOOD 
BAD CHARACTER RESULT ZERO 


C* ****** if * if # * * if * if if ** if * if if if if if if if * ** * * if if* ** * * * * ** ** **** ** ****************** * 


INTEGER INPUT! 1 ) f TYPE T TEMP , TEMP 1 , CHAR 
R E AL-4 ARGl 
I NT EGER *4 ARG2 
INTEGER .BLK t RP/PLtDP,E f CM 
LOGICAL SIN 

DATA MA5K1 tBLK f iP t RPf P.LtMI f OP^EtCM 
* / Z F0404040 ? 1 SM 1 * 1 ) 1 , • + ' * * 1 . 1 , •£% 1 , 1 / 

IEXP = 0 
IP = 1 


I SGN = 1 
ARGl = 0.0 


26 



ARG2 = 0 
SIN = .FALSE. 

SIGN = 1.0 
NUM = 1 
I START = K I R 
DO 20 i=I START, MAX 
KTR = I 

TEMP = INPUT (I) 

IF(TEMP.EO.BLK) GO TO 20 
GO TO (8,9,9,9,9,91 ,NUM 
8 IF(TEMP. NE.LP. AND. TEMP. NE.CM) GO TO 23 
NUM = 2 
GO TO 20 


9 I F ( TEMP . EQ.RP . OR . TEMP . EQ .CM ) GO TO 22 
GO TO (10, 10, 13, 15, 17, 19), NUM 

10 NUM = 3 

IF ( TEMP.NE.PL ) GO TO 11 
I F ( S I N .AND. TYPE.E0.2) GO TO 23 
SIN = .TRUE. 

GO TO 20 

11 IF(TEMP.NE.MI ) GO TO 13 

I F ( S I N .AND. TYPE. E 0.2) GO TO 23 
SIN = .TRUE. 

SIGN = -1.0 
GU TO 20 

13 IF(TEMP.EQ.E) GO TO 151 

I F ( TEMP . NE.DP ) GO TO 14 
I F ( TYPE .£0 .2 ) GO TO 23 
NUM = 4 

GO TO 20 

14 TEMPI = I SRL ( 24 , TEMP— MASK 1 ) 

IF( TEMPI .LT.O.OR.TEMP1.GT.9) GO TO 23 
I F ( TYPE . EO. 1 ) ARG 1 = ARG 1*1 0 . O+FLfJAT ( TEMP 1 ) 


IF! TYPE. EO. 2) ARG2 = ARG2*10+TEMP1 
GO TO 20 

15 IF ( TEMP. (ME. E) GO TO 16 
151 NUM = 5 

IF(TYPE.E0.2 ) GO TO 23 
GO TO 20 

16 TEMPI = I SRL ( 24, TEMP-MASK 1 ) 

I F ( TEMPI . LT .0 .OR, TEMPI .GT . 9 ) GO TO 23 
IF ( TYPE.EO. 1 ) ARG1 = ARGl+FLOAT ( TEMPI ) / 10 .0** I P 
IP = I P+1 
GO TO 20 

17 NUM = 6 

I F ( TEMP.NE.PL ) GO TO 18 
GO TO 20 

18 IF ( TEMP. NE, MI ) GO TO 19 
ISGN = -1 

GO TO 20 

19 TEMPI = I SRL ( 24 , TEMP— MASK 1 ) 

I F ( TEMPI .LT .0 .OR . TEMP l .GT ,9 ) GO TO 23 
I EXP = I 6XP*10+TEMP 1 
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20 CONTINUE 
ARG1 = 0.0 
ARG2 = 0 
RETURN 1 

22 I F ( TYPE • EQ • 1 ) ARC 1 = $ I GN*ARG1*10.0** ( I SGN* I EXP) 
IF{TYPE.E0.2 .AMD. STGN.LT.O.) ARG2 “ -ARG2 

IF ( TEMP.EO.CM ) RETURN 2 
RETURN 

23 ARG1 = 0.0 
ARG2 = 0 
CHAR = TEMP 
RETURN 3 
END 


C 

SUB ROUTINE C OM AND { I NP.UT ,N.M ,*,*,*,* ) 

IF ERASE, CALL WIPE 
IF DUMP, CALL DUMP 
RETURN 1 NAME NOT DEFINED 
RETURN 2 INVALID COMMAND 
RETURN 3 SAID STOP OR QUIT 
RETURN 4 ALL BLANKS FOUND 

C Jjf # ?;? # s|< j}c sj< sjc # # # i\< # ?|c ?',c sfc # >!' # # : ( 'c a>« # >\t # s}c 3$: jjs ;{c # # # # # # # % sjs ;|c jje # ;{i ;|c sje # sj: s|t sjc # ;',c ;',s # # sjc ;',c ;{c 5J 

DIMENSION INPUT { 190) ,NM{2 ) , I ERAS ( 2 ) 

DATA I L P , I R P , I E R A S , I ft L K , I QU I T , I GUM P / 1 { 1 ,« )* ,* 6« ,«RASE # , 

* 1 1 , ’QUIT* , 9 DUMP 8 / 

DATA I STOP/ 1 STOP 1 / 

C LOOK FOR PARENS 

DO 1 1=1,190 
KLP = I 

IF( INPUTU j.EO.ILP) GO TO 2 

1 CONTINUE 

C IS IT QUIT 

CALL NMP ACK ( I NPUT , 1 , 19 0 , NM , . TRUE . , £4, £4, £9 ) 

I F ( NM ( 1 ) .NE.IBLK) GO TO 4 

IF(NM(2) .NE.IOUIT .AND. NM (2) .NE . I STOP ) GO TO 7 
RETURN 3 

C IS IT DUMP 

7 IF { NM ( 2 } .ME. I DUMP) GO TO 4 
CALL DUMPY 
GO TO 8 

2 DO 3 I =KLP , 190 
KRP = I 

I F ( INPUT { I ) .EQ. IRP) GO TO 5 

3 CONTINUE 

4 RETURN 2 

C PACK ERASE 

5 INPUT! KLP) = I BLK 

CALL NMP ACK ( INPUT ,1 ,KLP , NM , . TRUE . , £4 , £4 , £9 ) 

I F ( NM ( 1 ) . NE . I ERAS ( 1 ) .OR. NM ( 2 ). NE . I ERAS ( 2 } ) GO TO 4 


28 



o o 


C COMMAND OK 

CALL NMPACK ( INPUT, KLP+1 ,KRP-1 ,NM, . FALSE . , €4, &4, £9 ) 
CALL F INONM ( NM , I NM , £6 ) 

CALL WIPE(INM) 

8 RETURN 

6 RETURN I 

9 RETURN 4 
END 


SUBROUTINE DUMPY 

TELLS USER WHAT IS IN BIG ARRAY 
COMMON/MATE/ NUM , I L AST , M ATR I X ( 1 500 ) 

IF(NUM.GT.O) GO TO 1 
PRINT 103 
GO TO 20 

1 I F ( NUM.GT . 1 ) GO TO 2 
PRINT 104 
GO TO 3 

2 PRINT 101, NUM 

3 J = 1 
DO 10 1=1, NUM 
NR = MATR I X ( J+3 ) 

NC = MATRIX (J+4) 

I F ( NR • EO . 1 .AND. NC.EO.l) GO TO 5 
C MATRIX PROPER 

PRINT 102, MATRIX! J ) , MATRIX! J+l ) ,NR,NC 
GO TO A 

C SCALAR 

5 PRINT 105, MATRIX! J) , MATRIX! J+l) 

A J = J + MATR I X ( J+2 ) 

10 CONTINUE 
20 RETURN 

101 FORMAT!' YOU HAVE DEFINED ',12,' MATRICES.') 

102 FORMAT (5X,2A4,5X, , HAS' , 13, * ROWS', 13,* COLUMNS.') 

103 FORMAT!' NO MATRICES DEFINED.') 

104 FORMAT!' YOU HAVE DEFINED ONE MATRIX.*) 

105 FORMAT ( 5X , 2 A4 , 5X , • IS A SCALAR.') 

END 
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c 


SUB R OUT I NE A R I TH ( K.K , I NM , IN Ml r I MM2 , N 1 , Ml , N2 , M2 , NC , MC , * ) 

c # $z >;< >\z >Js jjc >Js >;? 3|c ^ 5|c jJc 5|c ^ # s{c :{< # ;£ i{c # # # # i|c jJc.jJs # s|c :|t sjs.jjc ;|c jjc s}s jJc ;J c s|c sjc ;|c sj 


ARITH'S JOB IS TO DEFINE THE ANSWER MATRIX ENTRIES AND SIZE 
THIS ONE OSES FP NUMBERS IN THE COMMON 
KODE TABLE FOLLOWS 

1 — * 2 — + 3~ — 4-- -INV 5 — DET 6 — TRANS 7 — POWER 
8 — EIGEN 9 — DIAGADD 10 — D I AG SUB II — DI AGMULT 


s£ Jjc sjs ?;c ?;« # sjs sjc 5i< 5j< >;< sjs ;[c >Jc # i[c .#•# sjs ^ ;;; ;;c # # # 3j; %. »Js sfc ijc 5$: sjc ?;? sjc ^ #-3jc sjc sjc ?;< ^ # ;Jc ;Jc ^ >;< # ;|c >Jc 3|< 5|c t\< ?Jc ;[c ;Jc # j{c ?Jc # ;Js 

COMMON/MATE/ NUM , I L AST , RMATRX ( 1500) 

KODE = KK 
I N O = INM + 5 
I ND1 = INM I + 5 
NA = N 1 
MA ■= Ml 

GO TO { 1 3 * 1 3 f 1 3 * 4 ? 4 * 7 * 4 * 4 * 1 3 * 1 3 * 13 ) t KODF 
13 IND2 = I MM2 + 5 
NB = N2 


MB = M2 

I F ( MA • EO . 1 . AND • NA.EO. 1 ) GO TO 20 
I F ( MB . EQ • .1 •AND. NB.EO.l) GO TO 21 
I F ( KODE . L E . 8 ) GO TO 12 
PRINT 102 

102 FORMAT { 1 OPERATOR NOT APPLICABLE FOR MATRICES.') 

GO TO 11 

12 GO TO (1 *2 *3 *4*4*7 ? 4 * 4 ) , KODE 

1 C A L L M U L M A T ( R M ATR X ( IN D 1 ) 9 R M A T R X ( I M D 2 ) , R M A T R X ( I N D ) , N A , M A f N B , M B f & 1 1 ) 

NC = NA 

MC - MB 

GO TO 10 

2 CALL ADDMAT ( RMATRX < I N D 1 ) , RMATRX ( IND2 ) , RMATRX ( I NO ) , NA t M A t NB * M ft * £ 1 1 ) 

NC - NA 

MC * MA 

GO TO 10 

3 CALL SUBMAT ( RMATRX ( IND1 ) t RMATRX { IND2) , RMATRX { I NO ) , NA f MA , NB t MB v £ 1 1 ) 

NC = NA 

MC = MA 

GO TO 10 

4 IE(MA.EO.MA) GO TO 5 
41 PRINT 101 

101 FORMAT ( ' OPERATION NOT DEFINED FOR NON-SQUARE MATRIX.') 

GO TO 11 

5 I F ( KODF . EQ. 5 ) GO TO 6 
I F ( KODE • EO • 7 ) GO TO 8 
IF (KODE. EO • 8 ) GO TO 9 

51 C ALL MAT I NV ( RMATRX ( IND1 ) f RMATRX ( IMD) ,MA t £ll ) 

MC = MA 
NC = MA 
GO TO 10 

6 C ALL DETER ( RMATR X ( 1ND1 ) y RMATRX ( I ND ) t MA f DET *£11) 

RMATRX ( I ND) = DET 

MC .= 1 
NC = 1 
GO TO. 10 
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7 MC = NA 
NC = MA 

CALL TRANS! RMAT R X ( INDI ) ,RMATRX< IND) ,NA,MA) 

GO TO 10 

8 CALL EXPO! RMATRX! INDI ) , RMATR X ( I NO) , RMATRX ( 1400) f MA, INM2 ) 

MC = MA 

NC = MA 
GO TO 10 

9 IN 02 = IND + MA 

CALL E I GEN ! RMATR X ( 1400) , RMATRX ! I NTH ) , RMATR X ( IND) , RMATRX ! IND2 ) , MA ) 
MC = 2 
MC = NA 

10 RETURN 

11 RETURN 1 


SCALAR IS INVOLVED 
IT IS A 

20 I SCALE = 1 

nc - m 

MC = MB 
GO TO 22 
C IT IS B 

21 I SCALE = 2 
NC = NA 

MC ■= MA 

22 GO TO (23,24,25,4,4,7,4,4,26,26,26) , KOBE 

23 K0DE2 = 3 
GO TO 30 

24 KOBE 2 = 1 
GO TO 30 

25 KOBE 2 = 2 
GO TO 30 

C 

26 GO TO (27,28) * ISCALE 

27 I F (MB • ME • MB ) GO TO 41 
GO TO 29 

28 I F ( NA • ME ..MA ) GO TO 41 

29 KOBE 2 = KOBE - 5 

30 GO TO (32,33) , ISCALE 

32 CALL SCALAR (RMATRX! INDI ) , RMATRX! IND2 ) , RMATRX ! IND ) ,NB , MB , K0DF2 ) 
GO TO 10 

33 CALL SCALAR ( RMATRX (IMD2) , RMATRX! INDI) , RMATRX! IND) , MA , M A , K0DE2 ) 
GO TO 10 

END 


C 

SUBROUTINE AOOMAT ( A , B ,C ,NA , MA ,NB,MB ,* ) 
C PERFORMS MATRIX ADDITION 

LOGICAL PLUS 

DIMENSION A ( NA , MA ) , B { NB , MB ) , C ( NA , MA ) 
PLUS = •TRUE# 
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1 HI = NA 
Ml = MA 
HZ = NB 
M2 = MB 

I F ( Ml • EQ. M2 ) GO TO 2 
PRINT 101 

101 FORMAT { 8 MATRICES HAVE DIFFERENT NO. OF COLUMNS. OPERATION IGNORE 
*D. 1 ) 

RETURN 1 

2 IFCNl.EQ. HZ) GO TO 3 
PRINT 102 

102 FORMAT { 1 MATRICES HAVE DIFFERENT NO. OF ROWS. OPERATION IGNORED. 1 
* } 

RETURN 1 

3 IF (PLUS) GO TO 5 
DO 4 1=1 * Ml 

DO 4 J=1 ,N1 

4 C( Jt I ) = A(J,I) - B(JU) 

RETURN 

5 CONTINUE 

DO 6 1=1 9 Ml 
DO 6 J = 1 1 N 1 

6 C ( J , I ) = A ( J , I ) + B ( J * I ) 

RETURN 

ENTRY SUBMAT ( A* B , C * NA , MA , NB , MB , * ) 

C PERFORMS MATRIX SUBTRACTION 

PLUS = .FALSE. 

GO TO 1 
END 


SUBROUTINE MULMA T ( A » B , C * N A , MA , NB , MB ) 

PERFORMS MATRIX MULTIPLICATION 
PERFORMS A*B NOT REPEAT NOT B*A 
DIMENSION A<NA,MA) , B<NB r MB) ,C( : NAt MB) 

N 1 = NA 
Ml = MA 
N2 = NB 
M2 = MB 

I F ( Ml . EQ • N2 ) GO TO 1 
PRINT 101 

101 FORMAT ( 1 FIRST NO. COLUMNS NOT EQUAL SECOND NO. ROWS. 4 ,/, 
* 4 MATRIX MULTIPLY REQUEST IGNORED. 1 ) 

RETURN 1 

1 CONTINUE 

DO 2 1=1 1 M2 
DO 2 J = 1 * N1 

2 C ( J , I ) = 0 . 

DO 3 1=1 t Ml 
DO 3 J=1,M2 
DO 3 K = 1 f N 1 

3 C(K*J) = C(KtJ) + A<K,n*BUtJ> 

RETURN 

END 
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FORMS THE TRANSPOSE OF A 
SUBROUTINE TRANS ( A,C,NA,MA) 
DIMENSION A(NA,MA) ,C(MA,NA) 
DO 1 I = 1 » MA 
DO 1 0=1, NA 
1 C( I ,0 > = A ( 0 » I ) 

RETURN 

END 


PERFORMS MIXED ARITHMETIC 
A IS THE SCALAR 

KODE= 1*2,3 SCALAR ADD , SUB TRACT , MULT I PLY 
KQDE=4, 5 ,6 SCALAR DIAGONAL ADD , SUBTRACT , MULT I PLY 
SUBROUTINE SCAL AR < A , B ,C , NB , MB , KKK ) 

DIMENSION B ( NB ,MB ) ,C(NB,MB> 

KODE = KKK 
N = NB 
M = MB 

GO TO (5,10,15,20,20,20) , KODE 

5 DO 6 I - 1 , M 
DO 6 0 = 1 * N 

6 C( J, I ) = BIO, I ) + A 
GO TO 50 

10 DO 11 1 = 1, M 
DO 11 J = l, N 

11 C { J , I ) = B ( 0 * I ) - A 
GO TO 50 

15 DO 16 1 = 1, M 
DO 16 0=1, N 

16 C ( 0 , I ) = B ( 0 , I ) * A 
GO TO 50 

C DIAGONAL OPERATIONS 

20 DO 21 I = 1 , M 
DO 21 0=1, N 

21 C(0,I ) = B ( 0 , 1 ) 

KODE = KODE - 3 

GO TO (22,25,30) , KODE 

22 DO 23 I = 1 , M 

23 C ( 1,1) = C ( I , I ) + A 
GO TO 50 

25 DO 26 1=1, M 

26 C(I,I) = C(I,I) - A 
GO TO 50 

30 DO 31 I = 1 , M 

31 C ( I , I ) = C( 1,1) # A 
50 RETURN 

END 
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TAKES MATRIX A TO POWER IEXP, ANSWER IN B 
SUBROUTINE EXPOfA , B , C , N.A , I E ) 

DIMENSION A(NA,NA) ,B(NA,NA) ,C<NA,NA) 

IPOW = IE 
MA = NA 
DO 1 1=1, MA 
DO 1 J = 1 , MA 

1 B(j,n = 0*0 

I F ( IPOW.LE.O) GO TO 10 
I F ( I POW. EQ* 1 ) GO TO 8 
DO 2 I -1 , MA 
DO 2 J= 1 , MA 

2 C(J,I) = A ( J , I ) 

INITIALIZATION COMPLETE 
IPOW = IPOW - 1 
DO 7 1 1 = 1,1 POW 

PUT A*C IN B 
DO 3 1=1, MA 
DO 3 J = 1 , MA 
DO 3 K= 1 , MA 

3 B ( K , J ) = B ( K , J ) + A(K, I ) *C ( I , J » 

ARE WE DONE 

I F ( I I * EQ* I POW ) GO TO 20 

SAVE B IN C 
DO 4 1=1, MA 
DO 4 J= 1 , MA 

4 C(J,I) = B ( J , I ) 

NOW PUT 0 IN B 
DO 5 1=1 , MA 
DO 5 J=1,MA 

5 B (J ,1 ) = 0.0 

7 CONTINUE 
GO TO 20 

8 DO 9 1=1, MA 
DO 9 J = 1 , MA 

9 8 { J , I ) = A ( J , I ) 

GO TO 20 


10 m 11 I-1,MA 

11 B(T,I)«1.0 
20 RETURN 

MD 


INVERTS AA, INVERSE IN C 
SUBROUTINE MATINV ( AA,C ,MX , *) 
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C..... MATRIX INVERSION Hi GAUV- J0>GAM H f*pa f I-Vm 
LOGICAL MAT I N 

DOUBLE PRECISION A, ft, PftOu, Y 

01 ML NS I ON A A ( M X t M X ) * C ( MX , M / ) 9 / { |0, 10 ) 

DIMENSION Y< 10, 10) ,Z !2 , 10, 10 ) ,K 06) 

EQUIVALENCE ( Y t Z ) 

E Q U I V A L E NC E ( ABODE , EDCBA ) 

DATA DEL T t E PS / . 0 0 0 1 , . 1 E - 7 / 

MAT I N = .TRUE. 

GO TO 1 

C. ....ENTRY FOR OBTAINING DETERMINANT OF INPUT MATRIX 
ENTRY DETER! AA t C t MX,DETt» ) 

MAT I N= .FAL SE . 

C..... INITIALIZE ROUTINE AND TEST MX { = 0 ft 0 E ft OF MATRIX) 

1 M = MX 

ABODE = EDC BA 
DO 50 1 = 1 , M 
DO 50 J-l t M 
50 XU, I) = AA(J,I ) 

IFIM.GT.l) GO TO 5 

2 PROD=X (1,1) 

IF( PROD. NE. 0.0) GO TO 4 

3 I F (MAT I N ) GO TO 25 
DET * 0.0 

RETURN 

25 PRINT 101 

101 FORMAT ( 1 SINGULAR MATRIX. INVERSE NOT PERFORMED. 1 ) 
RETURN 1 

4 XM t l) = 1.0/ PROD 
GO TO 23 

5 PR0D=1.0 
MM=M-1 

DO 6 1=1, M 
M I > = I 
DO 6 J= 1 » M 

6 Y ( I , J ) = X ( I , J ) 

C.... .BEGIN BY FINDING LARGEST PIVOTAL ELEMENT 
DO 11 I =1 , M 
A = O.ODO 
DO 7 J = I , M 

IF ( DABS ( Y ( J , 1 ) ) .LE. A) GO TO 7 
A = DABS ( Y ( J , 1 ) ) 

L=J 

7 CONTINUE 

I F ( A .EQ. O.ODO) GO TO 3 
C..... REARRANGE ROWS AND ORDER ARRAY 
N=K ( L ) 

K(L)=K(I ) 

K ( I ) =N 
DO 8 J = 1 , M 
A = Y( I, J) 

Y ( I ,J )=Y(L,J ) 

8 Y { L , J ) = A 

C..... REDUCE PIVOTAL ROW 
A = Y ( 1,1) 

I F ( .NOT. MATIN) PROD=PROD*A 
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DO 9 J-l » MM 
9 Y ( I ,J >=Y< I * J + l ) /A 
Y!I,M) = 1 .000/ A 

C REDUCE REMAINING ROWS 

DO 11 L = 1,M 
IFd.EO.I ) GO TO 11 
A = Y ( L , 1 ) 

DO 10 N=1,MM 

Y ( L » N ) = Y ( L » N+1 ) -A*Y ( I ,N) 

10 IF(DABS< Y(L,N) ) ,LT, ( EPS*DABS ( Y ( L ,N+1 ) ) ) ) Y(L,N) =0.0D0 
Y(L,M)=-A*Y( I ,M) 

11 CONTINUE 

C UNSCRAMBLE INVERTED MATRIX 

DO 15 1=1, M 

I F ( K ( I ) . EQ . I ) GO TO 15 
PROD=-PROD 
DO 12 J = I , M 
I F ( K ( J ) . EO . I ) GO TO 13 

12 CONTINUE 
GO TO 3 

13 DO 14 L=1 , M 
A=Y(L, I ) 

Y ( L , I ) = Y ( L , J ) 

14 Y ( L , J ) =A 
K < J ) =K ( I ) 

15 CONTINUE 

IF( .NOT. MATIN) GO TO 23 

C OBTAIN ERROR MATRIX 

TEST-0.0 
DO 17 I = 1 , M 
DO 17 J=1,M 
R=0 .0 

DO 16 L = 1 , M 

16 R=R-Z ( 1,1, J ) *X ( I ,L ) 

IF( I.EQ.J) R = R + l.ODO 

TEST = AMAX1 (TEST,SNGL(DABS(R) ) ) 

17 Z ( 2 , I » J ) =R 
DO 19 1=1 ,M 
DO 19 J = 1 , M 
A = O.ODO 
DO 18 L= 1 , M 

18 A=A+Z ( 1 , I , L )*Z ( 2 ,L , J ) 

19 Z ( 1 , I , J ) =Z ( 1 , I , J )+A 

I F ( TEST .LE .DELT ) GO TO 21 
PRINT 102 

102 FORMAT!* MATRIX ILL-CONDITIONED. INVERSE MAY BE INACCURATE.*) 
C TRANSFER FINAL INVERSE 

21 DO 22 I = 1 » M 
DO 22 J = 1 , M 

22 C ( I , J ) = Z(1,I,J) 

23 IF { M. EO. 1 ) C( 1,1 ) = X €1,1) 

DET=PROD 

RETURN 

END 
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SUBROUTINE E I GEN ( A ,B ,.EVR , EV I , NN ) 

REAL*4 S 1 1 0 ) , MAX 
REAL*8 T 
COMPLEX R ( 10) 

DIMENSION A (NN,NN) ,B(NN,NN) , EVRI NN ) ,EVI (NN ) 
DATA EP/5.E-6/ 

M = NN 
NR = 1 
DO 1 1=1, M 
DO 1 J=1,M 
1 A I J , 1 ) = B I J , I ) 

SUM = 0. 

DO 210 1=1, M 
DO 210 K = 1,M 

210 SUM = SUM + ABS I A ( K , I ) ) 

EPS = EP*SUM/FLOAT(M*M> 

3 MM 1 = M - 1 
NM1 = M 
A N = NM1 

NM1 = N - 1 

5 IFIN.EQ.l ) GO TO 125 
MAX = 0. 

DO 10 I =1 » NMl 

I F ( MAX .GE. ABS I A IN, I ) ) ) GO TO 10 
MAX = ABS ( A ( Nt I ) ) 

IC = I 
10 CONTINUE 

IF(MAX.LE.EPS) GO TO 125 
IFUC.EQ.NM1) GO TO 40 
DO 20 K- 1 f N 
D = ALK,NM1) 

A ( K , NM1 ) = A(K, IC I 
20 A ( K 1 1 C ) = D 
DO 30 I = 1 » M 
D = A ( N M 1 , I ) 

A ( NM1 * I ) = At I C * I ) 

30 A ( IC » I ) = D 
40 D = 1 ./A (NtNMl ) 

DO 50 1=1, M 
SII) = A IN, I) 

50 AIN, I) = -AIN, I )*D 
AIN, NMl) = D 
DO 70 K=1 , M 
IFIK.60.NMl) GO TO 70 
E = A I N,K ) 

DO 60 1=1, NMl 

60 A I I ,K ) = A ( I ,K )+ A I I , NMl )#E 
70 CONTINUE 

DO 80 1=1, NMl 
80 All, NMl) = A I I , NMl ) *D 
DO 110 1=1 ,M 
T = 0.0 DO 
DO 100 K=1,NM1 
100 T = T + S(K)*A(K,I) 

110 A I NMl , I ) = T 

DO 120 I =NM1 ,MM1 

120 A I NMl , I ) = A I NMl , I ) + SII + 1) 

GO TO 4 

125 DO 130 I =N»M 
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130 S < I ) = -A ( N, I ) 

L = M - NM1 

CALL R00TXX<S(N),R(NR),L) 
IF(N.EQ.l) GO TO 150 
NR = NR + L 
M = NM1 
GO TO 3 

150 DO 160 1=1, NN 

EVR ( I ) = REAL (R( I ) ) 

EVI ( I ) = A I MAG { R ( I ) ) 

160 CONTINUE 
RETURN 
END 


SUBROUTINE TO SOLVE POLYNOMIALS WITH ROOTS ORDER. LE.2 
SUBROUTINE ROOTXX ( A * R,NN) 

COMPLEX CMPLX,C0NJG,CSQRT,A11,A12,A21,A22,A31,A32,S1,S2,X,DX, 

* TC,F0,F1 ,F2,R 

EQUIVALENCE I F0,A12) , ( F1,A22 ) , (F2,A31 ), ( TC,A11) 

* , (DX,A11) , (S2,A32) 

COMPLEX CZERO / (0.,0.) / 

COMPLEX CONE / { 1. ,0. ) / 

DIMENSION AM)fR(l) 

DATA MAX, EPS / 35,.5E-5 / 

N = NN 

IF(N.EQ.l) GO TO 150 
NP1 = N + 1 
DO 30 K-l , NN 
KK = K 

I F ( A ( N ) .NE.O.) GO TO 40 
N = N - 1 
30 R(K) = CZERO 
40 FN = N 

NM1 = N - 1 
FNM1 = NM1 
NP1 = N + 1 
KS = KK 

ST1 = A ( N-l ) / A( N ) 

ST2 = ST 1**2 - ( A ( N-2 ) *2 • ) / A ( N ) 

ASSIGN 110 TO Ml 
45 I F ( KK.GT.NN ) RETURN 
KM1 = KK - 1 
ASSIGN 70 TO M2 
C BEGIN ITERATION LOOP 

DO 120 L=1 ,MAX 
GO TO M2, (70,50) 
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50 A 1 2 = CONE 
A22 = CONE 
A32 = CONE 
DO 60 1=1, NM1 
All = A 1 2 

A12 = A ( I ) + X*A11 
A21 = A22 
A22 = A12 + X*A21 
A31 = A32 

60 A32 = A 22 + X*A31 
F2 = A31+A31 
FO = A(N) + X*A 12 

IF(REAL(FO>.EQ.O. .AND. A I MAG ( FO ) . EQ. 0 . ) GO TO 130 

51 = Fl/FO 

52 = S1*S1 - F2/F0 
GO TO 80 

70 X = CZERO 

51 = CMPLX ( ST 1 ,0 » ) 

52 = CMPLX ( ST2 » 0. ) 

ASSIGN 50 TO M2 

80 GO TO Ml, ( 110,90) 

90 DO 100 I =KS »KM1 
TC = 1.0/ { X-R ( I ) ) 

SI = SI - TC 
100 S2 = S2 - TC*TC 
110 TC = FN/S1 

DX = TC/ (1.0 + CS0RT(FNM1*(TC*S2/S1-1.0) ) ) 

X = X - DX 

IF(CABS(DX)/CABS(X) .LE. EPS) GO TO 130 
120 CONTINUE 
NN = KM 1 
PRINT 201 

201 FORMAT ( • **WARNING** EIGEN FAILED TO CONVERGE ON AN EIGENVALUE') 
130 ASSIGN 90 TO Ml 
R(KK) = X 
KK = KK + 1 

IF ( A I MAG < X ) . EQ.O. .OR. KK.GT.NN) GO TO 45 
R ( KK ) = CONJG(X) 

KK = KK + 1 
GO TO 45 

150 R ( 1 ) = CMPLX (-A(l ) ,0.) 

RETURN 

END 
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